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基因 序列 变异 与 遗传 分 化 
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摘要 : 采用 PCR 方法 测定 了 新 疆 东部 达 坂 城 (DB)、 吐鲁番 (TU)、 部 善 (SS)、 柳 树 泉 (HL) 和 哈密 东 郊 (HD) 
Huti (Phrynocephalus axillaris) 的 5 个 居 群 66 只 个 体 的 线粒体 ND4-tRNAIe 基因 序列 , 序列 比 对 得 到 838 
bp 位 点 ，29 个 变异 位 点 定义 了 17 种 单 倍 型 。 除 DB 居 群 外 ， 其 余 居 群 及 5 个 居 群 总 体 均 呈 现 单 倍 型 多 样 性 较 
f&1[(0.60050.113) A (0.922::0.012)],. 而 核 苷 酸 多 样 性 较 低 [(0.00082+0.00020) 夺 x 三 (0.00917+0.00037)]。 RRR 
分 析 支 持 5 个 居 群 形成 了 Clade I (DB+TU+SS) 和 Clade II (HL-HD) 两 个 分 支 。 单 倍 型 网 络 图 也 表明 ，5 个 居 群 
均 有 谱系 地 理 结构 。 分 子 变 异 分 析 (AMOVA) 表明 ，81.67% (P<0.01) 的 遗传 分 化 来 自 5 个 居 群 之 间 ， 而 距离 隔离 
分 析 (isolation by distance) 显示 各 居 群 之 间 的 遗传 分 化 (Fs 与 居 群 间 地 理 距 离 无 显著 关系 (r=0.5324，P=0.1196)。 
歧 点 分 布 (mismatch distribution) 分 析 和 Fu's Fs 检验 显示 ，TU 和 HL 居 群 可 能 经 历 近期 的 种 群 扩张 ， 而 SS 和 
HD 居 群 相对 稳定 。DB 居 群 特殊 的 栖息 生境 及 遗传 多 样 性 的 匮乏 揭示 该 种 群 可 能 处 于 瓶颈 期 。 推 测 更 新 世 以 来 ， 
这 些 区 域 的 地 质 、 气 候 变 动 引起 的 叶 城 沙 蜥 栖息 地 变化 及 遗传 漂 变 共同 导致 了 这 5 个 居 群 的 遗传 分 化 。 


关键 词 ， 叶 城 沙 蜥 ，ND4-tRNA'™， 遗 传 分 化 ， 新 疆 东部 
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Population Genetic Differentiation of Phrynocephalus axillaris in East 
of Xinjiang Uygur Autonomous Region Based on Sequence 
Variation of Mitochondrial ND4-tRNA'^" Gene 


LIJun"?, GUO Xian-Guang', WANG Yue-Zhao'' 


(1. Chengdu Institute of Biology, the Chinese Academy of Sciences, Chengdu 610041, China; 
2. Graduate University of the Chinese Academy of Sciences, Beijing | 100049, China) 


Abstract: A 838 bp fragment of mtDNA ND4-tRNA'™ gene was sequenced for 66 individuals from five populations 
(DB: Dabancheng, TU: Turpan, SS: Shanshan, HL: Liushuquan, HD: East district of Hami) of Phrynocephalus axillaris 
distributed in east of Xinjiang Uygur Autonomous Region. Seventeen haplotypes were identified from 29 nucleotide 
polymorphic sites in the aligned 838 bp sequence. Excluding DB, there were relatively high haplotype diversity 
[(0.60020.113) 5 1£(0.922:0.012)] and low nucleotide diversity[(0.000822-0.00020) X: x:(0.009172-0.00037)] in the 
whole and the rest four populations. The result of phylogenetic analysis showed that the haplotypes from five populations 
formed clade I (DB-TU- SS) and clade II (HL+HD). Network analysis further revealed the phylogeographical structure 
of these five populations. Analysis of molecular variance (AMOVA) indicated that the significant genetic divergence 
(81.6796, P«0.01) existed among populations, which was not obviously correlated with the geographic distances between 
populations inferred from isolation-with-distance analysis (7=0.5324, P=0.1196). The demography analysis implied that 
TU and HL might have experienced the recent population expansion; SS and HD might be relatively stable. With 
extremely low genetic diversity, DB might be influenced by the special environment and be experiencing the bottleneck 
effect. It followed that the significant genetic differentiation among the five populations of P. axillaris might be the 
consequence of habitat changes resulting from the geological events and climate oscillation since Pleistocene and genetic 
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生物 对 环境 与 气候 变化 的 响应 表现 为 表 型 、 生 
理 和 地 理 分 布 区 的 改变 以 及 适应 性 进化 或 者 新 的 
物种 产生 (Etterson & Shaw，2001)。 物 种 分 布 区 的 
变化 不 可 避免 地 影响 居 群 的 遗传 结构 (Hartl & 
Clark，1997)。 因 此 ， 对 居 群 遗传 结构 的 研究 可 以 
追溯 生物 居 群 的 历史 分 布 过 程 (Hewitt, 2000; Petit et 
al, 1997). 

vi ie — 2S UM EER CHR PCI 
WETS, AFIRE AI ro fit HUI BR 23-8 ETE 
旱 苑 漠 地 带 。 目 前 ， 基 于 沙 蜥 属 部 分 物种 的 居 群 遗 
传 结构 研究 发 现 ， 它 们 的 居 群 遗传 结构 明显 存在 坏 
境 变 迁 带 来 的 遗传 痕迹 (genetic signature)(Jin et al, 
2008; Jin & Liu, 2009; Liu etal, 2008; Wang et al, 
2009; Zhu, 2008) . I3 T PW (Phrynocephalus axillaris) 





































































































是 否 在 冰期 以 及 冰期 后 的 气候 变迁 中 存在 强烈 的 
种 群 动态 变化 ; (3) 天 山 山 脉 的 新 构造 运动 对 5 个 居 
群 遗传 分 化 的 影响 。 


1 材料 和 方法 


1.1 样品 采集 
叶 城 沙 晰 样品 来 自 新 疆 东 部 达 坂 城 、 吐 和 鲁 番 、 
部 善 、 柳 树 果 和 哈密 东 郊 5 个 居 群 , 共 66 只 个 体 (图 
1, X 1)。 取 新 鲜 肝 脏 或 肌肉 组 织 于 95% 酒 精 固定 ， 
标本 经 10% 福 尔 马 林 固定 后 转 70% 酒 精 保存 于 中 国 
科学 院 成 都 生物 研究 所 两 栖 爬 行动 物 标本 馆 。 
12 ”基因 组 DNA 提取 
取 适 量 已 固定 的 肌肉 或 肝脏 组 织 用 于 基因 组 
DNA 提取 。 采 用 常规 蛋白 酶 K 裂解 组 织 细胞 ， 酚 - 













































































































































































作为 我 国 特 有 的 仆 行 动物 ， 在 新 疆 境内 主要 分 布 于 
天 山 山脉 南部 的 荡 漠 、 半 亮 漠 地 带 (Zhao et al, 
1999)。 更 新 世 以 来 , 天 山 山脉 的 新 构造 运动 促进 了 
吐鲁番 -哈密 盆地 (简称 吐 哈 盆 地 ) 和 柴 窝 堡 盆地 之 
间 南 山 隆起 带 的 发 育 (Cao & He，1997) 及 天 山东 部 
支脉 (包括 吐 哈 盆地 中 的 了 壤 突 起) 的 隆 升 Jian， 
2001; Zhan，2002)， 对 吐 哈 贫 地 和 柴 帘 堡 贫 地 的 环 
境 产 生 了 很 大 影响 。 在 野外 考察 中 ， 我 们 发 现 新 疆 
东部 达 坂 城 地 区 的 叶 城 沙 蜥 体 色 、 背 部 斑纹 显著 不 
同 于 其 他 居 群 ， 而 且 砾 石 基底 的 栖息 生境 也 明显 有 
别 于 其 他 居 群 的 沙 质 、 沙 土质 基底 生境 。 这 种 形态 
分 化 是 否 具有 相应 的 谱系 遗传 分 化 是 令 人 感 兴趣 
的 问题 。 

线粒体 DNA(mitochondrial DNA, mtDNA) 以 较 
快 的 进化 速率 、 严 格 母系 遗传 、 几 乎 不 发 生 重 组 等 
特性 ， 能 够 较为 全 面 地 反映 居 群 内 和 居 群 间 的 遗传 
变异 (Avise et al，1987)， 因 而 被 广泛 用 于 群体 遗传 
学 和 分 子 系统 学 研究 。 在 本 研究 中 ， 我 们 拟 选 择 进 
化 速率 较 快 的 线粒体 ND4-tRNAIe 基因 为 分 子 标 
记 ， 通 过 其 序列 变异 对 新 疆 东 部 叶 城 沙 蜥 的 达 坂 城 
(DB)、 吐 鲁 番 (TU)、 部 善 (SS)、 哈 密 柳 树 泉 (HL) 和 哈 
密 东 郊 (简称 哈 东 , HD) 等 5 个 居 群 的 遗传 结构 进行 


















































































































































































































毛 法 (Sambrook et al，1989) 抽 提 ，95% 乙 醇 沉淀 获 
得 总 DNA。 提 取 完 毕 后 ， 用 50 pL 双 蒸 水 (ddH2O) 
溶解 总 DNA, 0.8% 琼 脂 糖 凝 股 电 泳 检测 后 用 于 PCR 
扩 增 。 
1.3 目的 片段 扩 增 及 序列 测定 

采用 Arévalo et al(1994) 中 的 引物 NDA 和 Leu 
进行 NDA-tRNA ^" 基因 片段 PCR. 扩 增 和 测序 ， 引 
物 由 北京 三 博 远志 生物 技术 有 限 责 任 公司 合成 。 
PCR 反应 体系 为 25 uL, 包括 10xbuffer( 含 MgCl>15 
mmol/L)2.5 uL, 2.5 mmol/L 的 dNTP 2 uL, 10 
pmol/uL 的 引物 各 1 uL, 5 U/uL 的 Taq 酶 0.2 uL, 
模板 DNA 1 uL(ZJ; DNAI0-— 100 ng). 最 后 加 灭 菌 
ddH;O 补 至 25 uL. PCR 反应 条 件 : 94C 预 变性 4 
min; 94C 变 性 40 s, 55C 退 火 40 s, 72C 延 伸 1 min, 
共 进 行 35 个 循环 ， 最 后 72C 延 伸 7 min; 反应 在 
PT200 PCR 仪 上 进行 。 PCR 产物 经 0.8% 琼 脂 糖 凝 
胶 电 泳 检测 后 送 上 海 英 骏 生物 科技 有 限 公 司 纯化 
测序 。 
14 数据 分 析 

采用 BioEdit 软件 (Hall，1999) 对 DNA 序列 进 
行 比 对 和 编辑 ， 同 时 辅 以 手工 矫正 。 用 MEGA4.0 
(Tamura et al，2007) 统 计 碱 基 组 成 、 变 异 位 点 数 等 
































































































































分 析 ， 探 讨 3 个 关键 问题 : (1)5 个 居 群 间 是 否 存 在 
明显 的 遗传 分 化 ， 具 有 不 同 的 单 倍 型 ，(2)5 个 居 和 群 























序列 特征 。 
根据 沙 蜥 属 的 系统 发 育 关 系 (Guo & Wang, 
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图 1 新 疆 东部 地 








区 叶 城 沙 蜥 采样 点 








Fig.1 The sampling localities of Phrynocephalus axillaris in east of Xinjiang Uygur Autonomous Region 


A: 南山 隆起 带 ，B: THRE. 


A and B represent the approximate site of Nanshan structural belt and Liaodun uplift, respectively. 


表 1 叶 城 沙 晰 采样 点 及 样品 信息 


Tab.1 Information of samples and their localities of Phrynocephalus axillaris 














种 群 采集 地 海拔 经 纬度 样本 量 标本 号 

Population Locality Altitude(m) ^ Coordinate Sample size Sample No. 
达 坂 城 N43.34205° 

DE Dabancheng 1198 E88.22084° 8 WGXG08527-534 
吐鲁番 N42.85885° 

TU Turpan —89 E89.190799 16 WGXG08209-213,WGXG08222-231,WGXG257 
a N42.84384° 

SS Shanshan 418 E90.19471? 11 WGXG08246-256 

HL BIROR 1085 BALA 15 WXJ09221-223,WXJ225-233,WXJ09236-238 
Liushuquan E92.90997° 72229, 7233, E 

HD 哈密 东部 828 e 16 WGXG08162-169,WGXG08178-182,WGXG08186-188 
East of Hami E93.79622° pod -le2, z 


2007ER KEW (Phrynocephalus mystaceus) 
(Genbank 登录 号 AY054055), "i Jibi PW Phrynoceph- 
alus helioscopus(GenBank 登录 号 AY054050) 和 变 
色 沙 蜥 (Pprymrocepjpnalus versicolor)(GenBank 登录 号 
AY054095) 为 外 类 和 群 ， 采 用 最 大 简约 法 (maximum 
parsimony，MP) 和 贝 叶 斯 推论 (Bayesian inference, 

BD 构建 单 倍 型 的 系统 发 育 树 。 在 PAUP*4.0b10 





















































分 支 节点 置信 和 度 (BP) 用 自 展 检验 (Bootstrap) 经 1000 
次 重复 抽样 得 到 。 用 jModelTest 0.1.1(Guindon & 
Gascuel, 2003; Posada，2008) 选 取 用 于 贝 叶 斯 分 析 
的 核 昔 酸 进 化 的 最 佳 模型 ， 在 MrBayes3.1(Ronquist 
& Huelsenbeck, 2003) 中 构建 贝 叶 斯 树 ， 以 后 验 概率 
(posterior probability，PP) 表 示 各 分 支 的 节点 置信 
度 。 起 始 树 为 随机 树 ， 以 3 条 热 链 、1 条 冷 链 的 马 

































































(Swofford，2003) 中 进行 最 大 简约 分 析 , 采用 启发 式 
搜索 (heuristic search), TBR(tree-bisection-reconnec- 


tion， 树 二 等 分 再 连接 ) 分 文 交换 搜索 最 大 简约 树 ， 





























尔 科 夫 蒙特 卡 罗 (Markov chain Monte Carlo)77 12:35 
ÍT 5x10f 4R, f$ 200 代 抽 取 一 棵 树 ， 最 初 1 万 棵 树 
作为 老化 样本 (burn-in) 被 舍弃 , 余下 的 1.5 万 棵 树 用 
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单 倍 型 网 络 图 。 
用 DnaSP 5.0 软件 (Librado 




















进一步 揭示 各 居 群 间 分 化 程度 。 











采用 Arlequin3.1(Excoffier et al; 
分 析 (analysis of molecular variance，AMOVA)， 探 
讨 叶 城 沙 蜥 的 遗传 结构 及 遗传 变异 的 分 布 情况 进 不 同居 群 中 的 分 布 见 表 2。 所 有 单 倍 型 中 ， 只 有 单 
而 作 靶 点 分 布 (mismatch distribution)(Rogers & 
Harpending，1992) 分 析 及 Fu's Fs(Fu, 1997)rp PERS 别 为 各 居 群 特有 








验 以 推测 种 群 是 否 经 历 扩张 。 
2 fh R 
2. 序列 变异 和 单 倍 型 分 布 





实验 获得 66 条 序列 (GenBank 登录 号 : 
GU952864-GU952921，GU952982-GU952989)， 比 X H11 和 Hl6。 


表 2 了 叶 城 沙 蜥 mtDNA ND4-tRNAr'™ 基因 单 倍 型 在 各 居 群 中 的 分 布 


于 计算 一 致 树 (consensus tree) 及 相关 参数 。 基 于 简 


约法 原理 用 TCS 1.21 软件 (Clement et al, 2000) 构 建 ND4 基因 的 3' 端 , 695—757 为 tRNA-His, 758—815 
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对 后 得 到 838 bp 的 片段 ， 其 中 1 一 694 为 mtDNA 





为 tRNA-Ser, 816—838 为 tRNA-Leu 的 5' 端 。ND4 


I 








& Rozas, 2009) 计 算 基因 的 终 1| 
各 居 群 的 单 倍 型 多 样 性 (站 、 核 苷 酸 多 样 性 m, 评价 BRET. A T C GF 
居 群 的 遗传 多 样 性 水 平 .Arlequin 3.1 软件 (Excoffier 
et al，2005) 中 计算 居 群 两 两 间 遗 传 分 化 指数 (FsD) 以 
用 IBD 1.52 软件 内 没有 发 现 操 
(Bohonak, 2002) 计 算 遗 传 分 化 是 否 遵循 距离 隔离 模 ， 括 19 个 简约 
式 (Isoation by distance, IBD), 与 地 理 距 离 呈 正 相 关 。 现在 密码 子 旬 
2005) 进 行 分 子 变异 tRNA 部 分 ( 占 13.8%); 密码 子 第 二 












































密码 子 为 TAA, 在 序列 内 部 未 发 现 终 止 
HEEN 35.296. 24.696. 
明显 的 反 G ffi, H. A+T 
含量 (59.8%) 明 显 高 于 C+G 含量 (40.2%)。 所 有 序列 
人 含 变异 位 点 29 个 ， 包 
i AVE SURI 10 个 单 突变 位 点 ， 分 别 出 











27.699. 12.699, WILH 
























































66 条 序列 共 定 义 了 17 个 单 倍 型 ， 各 单 倍 型 在 





倍 型 H4 为 SS f 
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第 三 位 ( 占 55.299) 
未 检测 到 变异 。 
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ITU 两 个 居 群 共享 其 余 单 倍 型 分 


ERE 1 个 单 倍 型 


Hl; TU 居 群 除 共享 单 倍 型 H4 外 ， 另 有 4 个 特有 





单 倍 型 ， 
































SS 居 群 除 












































有 4 个 特有 单 人 





H FRA H2 涵盖 了 该 居 群 的 10 个 个 体 ; 

ZW H4 外 , 也 有 3 个 特有 单 倍 型 ， 
其 中 单 倍 型 H8 出 现 频率 最 高 ; HL 和 HL 两 居 群 各 
所 含 个 体 数 最 多 的 单 倍 型 分 别 






































Tab.2 The distribution of mtDNAND4-tRNAIe" gene haplotypes in different populations of 
Phrynocephalus axillaris 





居 和 群 
Population(nl/n2/3) H1 H2 





H3 H4 H5 H6 H7 H8 H9 HIO 


单 倍 型 代码 Haplotype code 


H14 H15 HI6 H17 





DB(8/1/1) 7 
TU(16/5/4) 10 
SS(11/4/3) 

HL(15/4/4) 

HD(16/4/4) 

Total(66/17/15) 7 10 








The value in parentheses following popu 


22 ”系统 发 育 关 系 














依据 贝 叶 斯 信息 标准 (Bayesian information — 居 群 没有 形成 交互 单 




















3 3 1 1 3 5 1 4 


居 群 后 括号 中 数字 依次 代表 采样 数 、 单 倍 型 数 、 特 有 单 倍 型 数 。 





ation name represents number of sampling, haplo 


在 一 起 。 





clade I & ff; HL 





















































haplotypes, respectively. 


居 群 和 HD 居 群 ， 但 两 个 
^H HD 居 群 的 Hl4、H16 





criterion, BIC) 选 取 HKY+G 为 适用 于 贝 叶 斯 系统 发 。 RI H17 形成 一 个 亚 支 (PP=1.0，BP=95)。 


























育 推 论 的 核 音 酸 进化 最 佳 模型 。 

















简约 树 (B) 拓 扑 结构 大 体 一 致 (图 
居 群 形成 两 个 分 支 : clade I (PP=0.96，BP=65) 和 





























贝 叶 斯 树 (A) 和 最 大 

















从 单 倍 型 网 络 
2). "yb S (HD、HL) 的 单 倍 型 

















3) 可 知 ， 哈 密 地 区 两 居 群 
Ej TU. SS 以 及 DB 居 群 之 间 存 








在 较 多 未 检测 到 的 单 倍 型 ， 两 者 之 间 有 7 一 12 步 突 





cladeII(PP=-，BP=59)。cladeI 包 括 TU、DB 和 SS ” 4E, 与 图 2 中 的 cladeIl 和 clade 1 相对 应 ; TU JERE 


居 群 ， 其 中 DB 居 群 与 H4 之 外 的 TU 居 群 形成 一 

















文 且 互 为 姐妹 群 关系 , 共享 单 倍 型 H4 和 SS Ja FEX 











的 H2 与 相 邻 的 H3、H5、H6 仅 有 1 步 突变 ， 且 与 








DB 居 群 的 Hl 也 仅 有 1265828; 共享 单 倍 型 H4 与 
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H2 (TU) 
H3 (TU) 
H5 (TU) 
H6 (TU) 
H1 (DB) 
H7 (SS) 
Ho (SS) 
H4 (TU、 SS) — 
H8 (SS) — 
H17 (HD) 
H14 (HD) 
H16 (HD) 
上 H10 (HL) 
H12 (HL) 
H11 (HL) 
H15 (HD) 
H13 (HL) 


P. versicolor 






Clade I 





100 








0.62 
Clade II 








P. helioscapus 





P. mystaceus 
A B 


0.1 





图 2 ”基于 叶 城 沙 蜥 mtDNA ND4-tRNA" 单 倍 型 
构建 的 贝 叶 斯 树 (A) 和 最 大 简约 树 (B) 

Fig.2 Bayesian tree (A) and maximun parsimony tree (B) 
based on mtDNA ND4-tRNA'"" haplotypes of 
Phrynocephalus axillaris 

单 倍 型 代码 同 表 2. A 中 节点 数值 为 贝 叶 斯 后 验 概率 , B 中 为 1000 次 
复 抽样 得 到 的 自 举 值 ( 仅 显 示 三 50%)。 


Haplotype codes are the same indicated in Tab.2. Numbers at internodes are 






































Bayesian posterior probabilities for A and MP bootstrap values for B (only 
之 5$0% shown). 


SS 居 群 的 H7、H8 和 H9 分 别 有 1、2 和 2 步 突变 。 
此 外 ，HD 居 群 中 H15 并 未 与 该 居 群 其 他 单 倍 型 连 
在 一 起 ， 而 与 HL 居 群 中 H11 通过 1 步 突变 相连 ， 
类 似 于 HL 居 群 中 的 H10、H12 和 H13. 

2.3 ”遗传 多 样 性 和 居 群 遗传 结构 





























城 沙 蜥 5 个 居 群 中 除 DB 外 ， 其 余 四 个 居 群 和 
整个 研究 区 均 呈 现 单 倍 型 多 样 性 (有 ) 较 高 
[(0.600+0.113) ~ (0.92220.012)] ifti E EF B £ FE E (m) 
较 低 的 状 d  [(0.0008220.00020) ~ 
(0.00917+0.00037)]( 表 3). DB 居 群 因 只 含有 1 个 单 
倍 型 ， 因 此 ， 两 种 指数 值 在 该 居 群 中 均 趋 于 零 。 

AMOVA 分 析 显 示 ( 表 4, 将 叶 城 沙 蜥 5 个 居 群 
当 作 一 个 整体 时 ， 遗 传 变 异 主要 发 生 在 各 不 同居 群 
间 (81.67%，P<0.01)， 而 发 生 在 各 居 群 内 的 变异 仅 
占 18.33%(P<0.01); 其 他 分 组 分 析 中 尽管 得 到 所 占 
比率 较 高 的 组 间 变 异 ， 但 差异 显著 性 检验 均 表 明 其 
组 间 差 异 未 达到 显著 水 平 (P>0.05)。 居 群 间 遗传 分 
化 ( 表 5) 分 析 也 显示 : 两 两 居 群 间 的 遗传 分 化 均 达 到 
显著 性 水 平 (P<0.01)， 而 相对 较 小 的 遗传 分 化 分 别 
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图 3 ”基于 统计 简约 法 构建 的 叶 城 沙 蜥 单 倍 型 网 络 
Fig.3 The haplotype network Phrynocephalus axillaris 




















based on the statistical parsimony 
以 椭圆 表示 单 倍 型 , 其 代码 同 表 2。 各 单 倍 型 之 间 的 连 线 表示 一 步 变异 ， 
椭圆 大 小 约 表 示 该 单 倍 型 出 现 的 次 数 ， 小 空心 圆 代表 未 检测 到 的 单 倍 
型 。 


Each haplotype is represented by an oval and its code is the same indicated 










































































in Tab.2. Each mutation is shown as a line between each two haplotypes. 
The size of oval is proportional to corresponding haplotype frequencies. 


Small unfilled circles represent unsampled haplotypes. 


表 3 了 叶 城 沙 蜥 各 居 群 遗传 多 样 性 及 Fw's Fs 检验 
Tab.3 The genetic diversity and Fu's Fs test in different populations of Phrynocephalus axillaris 











居 群 单 倍 型 多 样 性 ATRE FENE Fu’s Fs 检验 显著 性 
Population Haplotype diversity (A) Nucleotide diversity (x) Fu's Fs test P-value 
DB a E Ex = 
TU 0.600+0.127 0.00188+0.00100 -0.29084 0.426 
SS 0.745+0.098 0.0021320.00029 0.50460 0.599 
HL 0.6002:0.113 0.000822-0.00020 -1.16127 0.088 
HD 0.7252-0.069 0.0041140.00078 3.15162 0.927 
Total 0.922+0.012 0.00917+0.00037 0.63457 0.629 
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表 4 时 城 沙 晰 mtDNAND4-tRNALe% 分 子 变 
Tab. 4 Hierarchical analysis of molecular variance (AMOVA) of mtDNA ND4-tRNA™ of Phrynocephalus axillaris 


变异 来 源 Source of Variation 


31 卷 


异 等 级 分 析 


变异 分 数 %variation © {Ë o value 显著 性 检验 P 





[DB, TU, SS][HL, HD] 


[DB,TU][SS][HL, HD] 


[DB,TU][SS][HL] [HD] 


[DB, TU, SS, HL, HD] 





Mw Me 
I I 


A tn 








H 
| 5i 





NS Ne 
I I 


居 和 群 后 








m 
zH 





Mw Me 
I I 


居 和 群 后 














hr fif a 























间 Among groups 


Among populations within groups 


EJ Within populations 


h] Among groups 


Among populations within groups 


F73 Within populations 


间 Among groups 


Among populations within groups 


居 群 内 Within populations 


Among populations within groups 


居 群 内 Within populations 


分 析 和 Fu's Fs 中 性 检验 结果 见 图 5 和 表 3. 明显 左 


de 5 叶 城 沙 蜥 各 居 群 间 遗 传 分 化 (Fsb) 及 差异 显著 性 


Tab.5 The pairwise genetic divergence Fst (blow 
diagonal) and its significance (above 
diagonal) between populations of 


Phrynocephalus axillaris 


TI HD Ji ff 


45.58 
39.17 
15.25 
64.23 
19.97 
15.81 
74.55 
8.14 
17.31 
81.67 
18.33 


Gct-0.45579 0.10564 
Gsc-0.84746 «0.01 
Gst-0.7197 «0.01 
Gct-0.64228 0.06733 
Gsc-0.84194 «0.01 
Gst-0.55815 «0.01 
Gct-0.7455 0.10119 
Gsc-0.8269 «0.01 
Gst-0.31983 «0.01 
Gsc-0.81671 «0.01 














T E] ELI 70 c 


























点 分 布 和 负 的 Fus Fs 值 显示 HL 和 




















TU 两 居 群 可 能 经 历 了 近期 的 种 群 扩张 ; 而 SS 居 群 
f 多 峰 的 歧 点 分 布 和 正 的 Fu’s Fs 值 则 说 明 





























种 群 处 于 相对 稳定 的 状态 。 





3.4. 居 和 群 遗传 结构 和 遗传 多 样 性 























DB TU SS HL HD 
DB P<0.01 P<0.01 P<0.01 P<0.01 3 i it 
TU 0.43983 P«0.01 P«0.01 P«0.01 
SS 0.87722 0.81655 P«0.01 P«0.01 系统 发 A 
HL 0.94658 0.88286 0.89398 P<0.01 
HD 0.79490 0.79555 0.79400 0.58881 1.00 

. 080 . 

出 现在 DB 与 TU 居 群 (Fsr=0.43983) 和 HD 与 HL 居 woof e 
群 (Fsf=0.58881) 之 间 ， 其 余 各 居 群 之 间 的 遗传 分 化 z ; 























的 关系 (=0.5324，P=0.1196)。 
2.4 种 群 历史 动态 
叶 城 沙 蜥 TU、SS、HD 和 HL 居 群 的 歧 点 分 布 


频率 Frequency 





柳树 泉 
"Hi 





指数 相对 较 高 ; 而 最 大 的 遗传 分 化 出 现在 HL 和 DB 
居 群 之 间 。 此 外 ，IBD 分 析 表 明 ( 
出 现 的 遗传 分 化 (Fs 与 





图 和， 各 居 群 之 间 














居 群 间 地 理 距 离 没 有 显著 




















xo …。-… 观 测 值 Obs 
居 群 期 望 值 Exp 
EN XI Obs 
期 望 值 Exp 

10 15 20 

















分 析 支 持 新 疆 东部 叶 城 沙 蜥 5 个 居 
r=0.5324 
e P-0.1196 





100 200 300 400 500 
地 理 距离 Geographic distance (km) 


图 4 叶 城 沙 蜥 距离 隔离 模式 分 析 


Fig.4 Isolation-by-distance analysis of Phrynocephalus 


axillaris 


0.4 


配对 差异 Pairwise differences 


EEERCUU 观测 值 Obs 


一 一 一 期 望 值 Exp 
T sd T 
"m n 





图 5 叶 城 沙 晰 部 分 居 群 单 倍 型 歧 点 分 布 分 析 
Fig. 5 The mismatch distributions of the haplotypes of Phrynocephalus axillaris in SS, TU, HD and HL 
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4 其 E 俊 等 : 


à 





群 形成 clade I (DB+TU+SS) 和 clade II(HL+HD) 两 个 








分 支 (图 2)。 同样 , 单 倍 








型 网 络 





图 也 文 持 了 这 种 遗传 





结构 : DB 居 群 与 TU 居 群 中 除 

















H4 外 的 其 余 单 倍 型 














关系 较 近 , 而 后 者 与 SS 居 和 群 








D E 
EX 


EEH H4 暗示 TU 和 
































S 


un 





居 群 个 体 之 间 存 在 一 定 程 








度 的 扩散 ;HL 居 群 和 
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TU 和 HL 两 个 居 群 在 其 演化 历史 上 可 能 经 历 过 种 
群 扩张 。 结 合 前 人 的 研究 (Grant & Bowen, 1998; 
Li etal, 2003; Liu etal, 2008; Xu etal, 2005), 
我 们 推测 居 群 内 出 现 单 倍 型 多 样 性 (h) 较 高 , 而 低 核 
音 酸 多 样 性 四 相 对 较 低 可 能 是 由 于 居 群 在 较 近 时 






































TER 








HD 居 群 (H15 除外 ) 之 间 相 对 较 少 的 缺失 单 倍 型 ( 仅 
4 个) 表明 它们 相对 较 近 的 关系 。AMOVA 分 析 发 现 











高 比例 的 遗传 变异 出 现 





n 


在 不 同 








居 群 间 (81.67%)， 表 





明 各 居 群 之 间 具 有 显著 的 遗传 分 化 。 居 群 两 两 间 的 

















Fst 值 也 再 次 说 明了 各 居 群 
管 DB 与 TU 和 HL 与 HD 





























昌 上 其 体 的 分 化 程度 。 尽 
大群 具有 相对 较 小 的 Fst 

















fü, 但 这 种 差异 仍然 显著 (P<0.01)。DB 居 群 的 个 体 


















































体 色 、 背 部 斑纹 及 栖息 生境 特征 不 同 于 其 他 四 居 
群 ， 而 与 四 居 群 间 显著 的 Fst 值 也 说 明 该 居 群 在 遗 





传 上 明显 有 别 于 其 他 居 群 。 基 于 居 群 间 遗 传 分 化 程 
度 和 采样 居 群 近似 线 状 的 分 布 所 做 的 距离 隔离 分 
析 (BD) 显 示 各 居 群 之 间 的 遗传 分 化 (fs) 与 居 群 间 











地 理 距离 没有 显著 的 关 
































2o 





建 群 种 越 小 越 容易 受到 遗传 漂 变 的 影响 ， 其 主 


要 结果 是 引起 遗传 多 样 性 的 丧失 (Stone，2002);， 而 























栖息 地 片段 化 所 导致 的 基因 
传 多 样 性 的 降低 (Xu et al，2005)， 增 加 居 群 间 的 遗 
传 分 化 (Frankham，1995)。5 个 居 群 之 间 出 现 明显 











ji Et IRE Se SOR 























间 范 围 内 经 历 的 瓶颈 后 的 快速 扩张 ， 有 利于 促进 单 
倍 型 多 样 性 的 增加 ， 但 还 没有 足够 的 时 间 积 累 其 核 
ERZ EEE. 
3.33 5 个 居 群 显著 遗传 分 化 的 成 因 

本 研究 中 , 5 个 居 群 分 别 位 于 叶 重 番 - 哈 密 盆地 
(简称 吐 哈 倪 地 ) 及 其 西端 的 柴 窜 堡 倪 地 达 坂 城 地 
区 。 吐 哈 盆 地 作为 天 山 山 脉 内 部 最 大 的 山 间 前 陆 贫 
地 ， 自 西向 东 分 别 为 吐鲁番 凹陷 、 了 墩 突起 和 哈密 
I"IA(Zhao et al, 1992); 有 旦 以 南山 隆起 带 与 柴 窜 堡 
盆地 相隔 (Allen etal, 1993). HD 和 HL 居 群 位 于 哈 
ak", TU 和 SS gv Tu RSS, DB 居 群 
ALT 3E 83 £6 dà Hh Zt ARIRE 1), 各 居 群 栖 
息 生 境 存 在 较 大 差别 。 更 新 世 以 来 。 天 山 山脉 新 构 
造 运动 促进 了 南山 隆起 带 的 发 育 (Cao & He, 1997) 
和 天 山东 部 支脉 中 了 墩 突起 的 隆 升 Jian，2001; 
Zhan，2002)。 南 山 隆 起 带 的 发 育 可 能 构成 DB 居 群 












































































































































同 TU 居 群 交流 的 障碍 ; 而 对 DB 居 群 而 言 ， 由 于 




















的 遗传 分 化 ， 但 系统 树 较 浅 的 分 文 结构 、 网 络 图 的 





连续 以 及 各 居 群 ( 除 DB 外 ) 较 高 的 单 倍 型 多 样 性 (及) 





和 较 低 的 核 苷 酸 多 样 性 (四 说 明 : 不 同居 群 很 有 可 能 























在 演化 过 程 中 受到 黄 基 者 效应 和 遗传 漂 变 的 影响 ， 

形成 了 遗传 多 样 性 低 的 现状 ， 而 在 适应 各 自 栖息 生 

境 过 程 中 逐渐 形成 了 居 群 的 特有 单 倍 型 。 

3.2 ”瓶颈 效应 和 近期 扩张 
DB 居 群 8 个 个 体 均 为 同一 单 倍 型 HL， 该 点 位 

于 天 山南 山 隆 起 带 西 侧 较 低 处 ， 海 拔 相 对 于 其 他 居 

群 最 高 (为 1108 m)， 基 底 表 层 满 布 砾石 ， 且 距 盐湖 
















































































不 足 1 km, 这 些 生境 特征 





他 居 群 的 沙 质 和 沙土 质 














明显 有 别 于 隆起 带 以 东 其 
生境 。 因 此 ， 推 测 DB Ju f 























单 倍 型 多 样 性 的 匮乏 是 











应 的 影响 。 单 倍 型 网 络 
中 心 的 辐射 状 结构 ， 并 








有 和 





于 受到 特殊 生境 瓶 须 效 
图 中 TU 居 群 形成 以 H2 为 

















SS 居 群 共享 单 倍 型 

















H4; HL 居 群 类 似 的 网 络 结构 以 及 中 央 单 倍 型 H11 


与 HD 居 群 的 H15 





仅 以 一 








步 突变 相连 ， 表 明 














两 居 群 可 能 发 生 过 扩张 事件 。 歧 点 分 布 分 析 中 显著 
的 单 峰 (图 5) 和 负 的 Fu's Fs 值 ( 表 3) 也 支持 






































自身 缺乏 有 效 种 群 的 调节 ， 从 而 容易 受到 黄 基 者 效 
应 和 遗传 漂 变 的 影响 ， 致 使 现 生 居 群 遗传 多 样 性 极 
低 。 了 墩 突起 可 能 阻隔 HL 居 群 和 SS 居 群 间 的 扩 
散 ， 使 其 形成 各 自 特有 单 倍 型 。 

第 四 纪 的 冰期 、 间 冰期 回旋 对 北半球 很 多 物种 
的 分 布 格局 的 形成 都 有 着 显著 的 影响 (Hewitt， 
1996, 1999); 冰期 的 干 冷 , 促进 了 欧 漠 及 沙化 环境 
的 扩展 (Wang，1997; Wu & Guo，2000; Wu et al, 
2002), 而 间 冰 期 相对 暖 湿 的 气候 有 利于 促进 绿洲 发 
育 ， 绿 洲 的 形成 可 能 造成 适宜 沙 蜥 生存 的 误 漠 环境 
的 片段 化 ， 进 而 对 5 个 居 群 的 遗传 分 化 产生 影响 。 
Liu et al(2008) 在 对 若尔盖 湿地 青海 沙 蜥 红 原 亚 种 
(Phrynocephalus vlangalii hongyuanensis) 的 研究 中 
岂 提 及 间 冰 期 暧 湿 的 气候 使 得 邻近 高 山上 的 雪 水 
融化 进入 低地 而 形成 的 沼泽 、 草 甸 阻 隔 了 不 同 沙丘 
间 青 海 沙 蜥 的 基因 交流 ， 最 终 使 得 不 同 地 理 单元 间 
出 现 显 著 的 遗传 分 化 。 因此 , 我 们 推测 更 新 世 以 来 ， 
南山 隆起 带 、 了 墩 突 起 以 及 各 居 和 群 栖息 地 之 间 绿 洲 
的 演化 使 得 适合 叶 城 沙 蜥 生存 的 充 漠 环境 片段 化 ， 
进而 导致 5 个 居 群 呈现 出 显著 的 遗传 分 化 。 
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